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An asymptotic restoration procedure is applied for ana- 
lyzing bound-state overlap functions, separation energies and 
single-nucleon spectroscopic factors by means of a model one- 
body density matrix emerging from the Jastrow correlation 
method in its lowest order approximation for ^^'O and ^"Ca 
nuclei . Comparison is made with available experimental data 
and mean-field and natural orbital representation results. 

PACS number(s): 21.10.Gv,21.60.Cs,21.60.Jz 



Both the single-nucleon spectroscopic factors and the 
overlap functions have attracted much attention in ques- 
tions concerning the interpretation of recent (e, e/p), 

{d,^ He) and (7,p) experimental data (e.g. [|l|-fi5[. The 
growing interest is motivated in principle by the possibil- 
ity to clarify the limitation of the nuclear mean-field pic- 
ture. For instance, the relatively low values of the spec- 
troscopic factors deduced from these experiments show 
clearly the importance of the short-range correlation ef- 
fects in nuclei and the necessity of detailed investigations 
of the high-momentum components of the nucleon wave 
function which cannot be included within the mean-field 
approximation |l^-|20||. 

The underlying relationship between the differential 
cross-section and the structure of the nuclear wave func- 
tion is often empirically analyzed using the plane-wave 
impulse approximation (PWIA). For instance, in this ap- 
proximation the (e, e/p)-reaction cross-section for a tran- 
sition to a specific state with quantum numbers a in the 
residual nucleus has the following form (see e.g. 



(e,e/p) 



d^lsQfd^pdEp 
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The first term /C is kinematical factor , tJep is the off- 
shell electron-proton scattering cross-section |^ and the 
nuclear structure component | 0ct(k) p is the squared 
Fourier transform of the overlap function between the 
ground state of the target nucleus "^^^^ and the final 
state of the residual nucleus 0,^,|T6|: 



Mr) = {^i""-'^ I a(r) | ^^^'^) 



(2) 



a(r) being an annihilation operator for a nucleon with 
spatial coordinate r (spin and isospin coordinates are 
not put in evidence). The overlap functions (^) are not 
orthonormalized. Their norm defines the spectroscopic 
factor of the level a 



Sa = (0Q I (t>a) 

and the normalized overlap function 
Mr) = S~'/^ Mr) 
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Usually, 4'a{'r) is calculated from an empirical Saxon- 
Woods potential with a distinct potential radius for each 
separate transition a. Quantitative estimates are then 
deduced by fitting both the potential radius and the 
spectroscopic factor Sa in order to obtain a good agree- 
ment between the experimentally measured cross sections 
ida/dri)^^p and those predicted from appropriate calcu- 
lations for the reaction process. 

The full theoretical description of the experiments 
mentioned above has many components. We should like 
to mention among them the proper account of the reac- 
tion mechanism, of the distortion effects (including the 
distortion due to the final state interaction), of the me- 
son exchange currents contributions, the study of the 
A-dependence and others (see e.g. |^,^). Obviously, 
however, the theoretical estimates of the overlap func- 
tions ^ and the spectroscopic factors (||) are of crucial 
importance for the adequate evaluation of recent (e, efp), 
{d,^ He) and (7,^5) experiments. The problem is that 
the normalized overlap functions (^ cannot be identi- 
fied with a phenomenological shell-model single-particle 
wave function especially for energies farther from the 
Fermi energy, sometimes even within the valence shell 
|]l6| . It is known that generally, the independent -particle 
shell model cannot explain the fragmentation or spread- 
ing of the hole strength |p6| , p2| , po| . This is because, due 
to the residual interaction, the hole state of the target 
nucleus is not an eigenstate of the {A — l)-nucleon sys- 
tem and its strength is distributed over several states of 
the final system. Possible modifications going beyond the 
uncorrelated shell-model approximation quickly become 
rather involved. 

Recently, it has been shown that the knowledge of 
the ground-state one-body density matrix of the target 
nucleus is sufficient to determine, at least in principle, 
the overlap functions, spectroscopic factors and separa- 
tion energies of the bound {A — l)-particle eigenstates. 
The aim of the present paper is to apply the procedure 
suggested in [ p7| using the model one-body density ma- 



trix |28|-|30|| in which the short-range correlation terms of 
the Jastrow correlation method are taken into account. 
The resulting quantitative estimates allow us to make 
instructive conclusions for the properties of the overlap 
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functions in comparison with the associated shell-model 
orbitals and the natural orbitalsj 3l[ which are of frequent 
interest in this context ||l^,|l^,|0 1 . 

The exact one-body density matrix associated with the 
ground state vjf^"^) of the target nucleus with A nucleons 
is defined as 



p(r,r/) = (*(^) I at(r) a(r/) | *(^)) 



(5) 



Inserting the complete set of eigenstates ^i^^ '^'^ of the 
residual {A — l)-nucleus, eq. (^) reads 

p(r,r/) = J2 '^"(^) '^"(^') ^*c.ir) Mr') , (6) 



where (f>a and (pa are the overlap functions (0) and (^), 
respectively, Sa is the spectroscopic factor and the 
summation implicitly includes also the continuum states 
associated with all scattering channels of the {A~l) sys- 
tem. 

The one-body density matrix has a quite similar form 
in its natural orbital representation |3l| 



(7) 



where the natural orbitals 4'a are defined as a com- 
plete orthonormal set of functions which diagonalize the 
one-body density matrix (j^) with eigenvalues called 
natural occupation numbers. Properties of the over- 
lap functions and the natural orbitals are considered 
for various many-body systems, such as atomic nuclei 
i^,p8|j29p|,|20|,|3^P] and liquid drops of ^He^. 



In the case of a spherical ground state \I>(^) with 0+ 
angular momentum and parity, each eigenstate ^i"^ is 
characterized by the 'single-particle' quantum numbers 
Ijm, i.e., a = nljm with n being the number of the Ijm- 
state. The overlap functions and the natural orbitals 
then factorize into radial and spin-angular parts 



(flair) = (pnl]{r) Fy„i(rj,cr), 

il^air) = i^nij{r) Yijm{^,a-), 



(8) 



where Yij^i^, ct) = x X^^^i'^)]m ^ is the spin 

variable. Using eq. (0), the one-body density matrix 
reads 

p(r,r/)=^ py(r,r/) ^ Y*^{^,a) Yi,„,{n/,a) . (9) 

Ij ma 

Due to the spherical symmetry Snij , N^ij and the radial 
contributions pij (r, rf) entering the one-body density ma- 
trix (|^) do not depend on the magnetic quantum number 
TO. From eqs. and (|^) then follows that in each Ij- 
subspace the spectroscopic factor Snij is smaller than the 
largest natural occupation number N^^j^ with the same 
Ij, i.e., 



Snl, < N, 



nlj 



(10) 



The procedure of is based on the generally ac- 
cepted asymptotic behavior of the neutron overlap func- 
tions associated with the bound states of the {A — 1) 
system 



i'nioir) Cnij exp(-fc„;j- r) / r , 



where 



AA-l) 



- E, 



(11) 



(12) 



{A-l) 
nlj 



E, 



(A) 



depends on the separation energy = E, 
Obviously, the higher excited states have faster decay. 
Therefore, for large values rf — a —^ oo, eqs. (|) and (|) 
lead to the asymptotic relation 

pij{r,a) (f'noiji'r) x Cnoij exp(-fc„„(ja)/a , (13) 

where (pnoijif) is the radial part of the lowest overlap 
function in the Ij subspace considered. The unknown 
constant Cngij can be derived from the asymptotic be- 
havior of the diagonal part pij (r, r) since 



Cn„ij P exp(-2fc„o;ja)/a^ 



(14) 



By means of eqs. (|T^) and ( [T^ ) one can derive the low- 
est (bound state) overlap function with quantum num- 
bers Ij and radial part 



Miji^) 



Pijir,a) 



Cnaij exp{~k„gija)/a 
as well as the associated separation energy 



(15) 



enoij = kl^ij/^m (16) 



and the spectroscopic factor 



Snolj — 



Mlj) ■ 



(17) 



One can repeat the above procedure for the second bound 
state with the same multipolarity (if it exists) after sub- 
tracting the contribution of the lowest eigenstate. The 
result is 



Mi]{r) 



pij{r,a) - (f)„oij{r) Mijia) 
Cmi] cxp(-fc„iya)/a 



(18) 



and expressions similar to (|T^ and ( [TtI ) for the sepa- 
ration energy tmij and the spectroscopic factor Smij, 
respectively. The restoration procedure can be contin- 
ued and one is able to analyze all bound states of the 
{A — l)-particle system once the one-body density ma- 
trix of the A-particle ground state is known. In the case 
of proton bound-states some modifications due to the 
Coulomb asymptotic behavior of the overlap functions 
have to be taken into account. 

The present calculations of the bound-state overlap 
functions, separation energies and spectroscopic factors 
have been performed applying the recipe (O) - (18) to a 
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model one-body matrix ps]-^ obtained within the Jas- 
trow correlation method in its low-order approximation 
for ^^O and ''°Ca nuclei . The model is based on har- 
monic oscillator single Slater determinant and Gaussian- 
like state-independent correlation factor. Although the 
resulting density matrix has a simple analytical form, it 
is physically significant that the short-range correlations 
are incorporated in it to a large extent. In addition, its 
natural orbital representation is well investigated [p8|-p0[ . 

Here we should like to mention that the procedure sug- 
gested in requires accurate values of the one-body 
density matrix at large r. In principle this limits the 
practical application of the method. In our opinion, how- 
ever, the analytical expressions of the one-body density 
matrix obtained within some correlation methods give a 
basis for an easier numerical application of the proce- 
dure. This is the case in the present work. Using the 
model one-body density matrix from psj-^of we have to 
avoid, however, another difficulty arising from its Gaus- 
sian asymptotic behavior. This is achieved by applying 
the recipe (^3|) - (|l^) not for a single asymptotic point a 
but within an asymptotic region {aL,au) bracketing the 
point a and sustained after the point where the diagonal 
element pij{r,r) is less than 10 percent from its maxi- 
mal value. We are looking for such a radial contribution 
Pnoij{r, rf) = (/Jnouir) (j^noijirf) whose diagonal part is less 
than or equal to pij (r, r) at each point r and which mini- 
mizes the trace Tr [{pij — Pnoij)'^] ■ Then a^, a and ajj as 
well as the unknown Cngij and kn^ij are uniquely deter- 
mined by the requirement that the overlap function (|l5| ) 
satisfies eqs. (^3|) and (^4|) simultaneously with minimal 
least-squared deviation within the region {aL,aij). We 
should mention that the procedure suggested is not the 
unique one, but this problem does not exist when a real- 
istic one-body density matrix with a correct exponential 
asymptotic behavior is considered. 

We have performed the above numerical procedure sep- 
arately for each set of quantum numbers nl (the model 
does not split the states with respect to j — I ± 1/2). It 
leads to predictions for the neutron separation energies 
Eq., spectroscopic factors Sa and the overlap functions (j>a 
which are given in Table I and Fig. 1, respectively. 

It is seen from Table I that the calculated separation 
energies are in acceptable agreement with the self- 
consistent Hartree-Fock (HF) results and the avail- 
able experimental data. The calculated spectroscopic 
factors Sa, however, differ significantly from their mean- 
field values. Due to the short-range correlations a de- 
pletion of the states below and a filling of the states 
above the Fermi level results. At the same time the cal- 
culated values of Sa are consistent with experimentally 
deduced spectroscopic factors |6|. In general, the values 
of Sa emerging from the present restoration procedure 
are larger than the experimental ones. This fact is most 
probably related to the crude approximation for the den- 
sity matrix used. The reason is the same for the larger 
value of S2SI/2 in comparison with the spectroscopic fac- 



tor of the lower Is-state in *°Ca. 

In Table I we compare Sa also with the natural occu- 
pation numbers Na derived after diagonalyzing the same 
model one-body density matrix |2^-|3C|] . The comparison 
shows that our numerical procedure satisfies the general 
requirement ([lo|). The trend of the calculated spectro- 
scopic factors Sa follows that one of the natural occu- 
pation numbers. This result becomes more transparent 
realizing that the overlap functions 0q are rather close 
to the natural orbitals tpa as it is seen from Fig. 1. 

From Fig. 1 it can be also seen that all three func- 
tions, the overlap, mean-field and natural orbital wave 
functions, are quite similar for the hole states in nuclei. 
This justifies the use of shell-model orbitals instead of 
overlap functions within PWIA calculations (|^) for such 
kind of nuclear states. This approximation, however, is 
no longer valid for the particle nuclear states where the 
mean-field wave-functions significantly differ from the 
overlap functions (see the 1 /-state in Fig. 1). The lat- 
ter take some intermediate position between the natu- 
ral orbitals and the HF wave-functions. It should be 
stressed that our model one-body density matrix is com- 
pletely different from the Hartree-Fock (HF) one. It has 
been demonstrated in |^,^,^ that due to the short- 
range correlations the correlated particle-state orbitals 
are much more localized than the particle-state mean- 
field single-particle wave functions. This is the reason 
why the HF orbitals associated with the particle-states 
go further out than the overlap functions (or the natural 
orbitals) associated with the correlated one-body den- 
sity matrix. The place of the correlated particle-state 
asymptotic region is affected by the SRC while the HF 
particle-state orbitals have a larger spread although the 
orbit is more strongly bound. 

The instructive conclusion is that neither natural or- 
bitals nor shell-model wave-function can be used as 
particle-state overlap functions in the theoretical anal- 
ysis of the experimental data. The present restoration 
procedure gives a possible solution of the problem if it is 
applied to some realistic ground-state one-body density 
matrices. 

Concluding, we have demonstrated in this paper the 
possibility for restoring the separation energies, spectro- 
scopic factors and overlap functions for bound {A — 1)- 
particle eigenstates on the basis of the ground-state one- 
body density matrix of the target j4-particle system. Al- 
though we have used quite crude approximation for the 
one-body density matrix [^8|-^ , the asymptotic restor- 
ing procedure leads to acceptable quantitative re- 
sults. Thus one obtains a method for estimating such 
important quantities as spectroscopic factors and over- 
lap functions which is supplemental to the more involved 
approaches |]l6|. For this purpose, one has simply to 
apply the present restoring procedure to more sophisti- 
cated one-body density matrices as for example the ones 
emerging from Brueckner-Hartree-Fock varia- 
tional Monte-Carlo or simplectic model calcula- 
tions The resulting bound-state spectroscopic fac- 
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tors and overlap functions will have more realistic prop- 
erties and can be used for reliable description of the char- 
acteristics of the (e, eip), {d,^ He) , (7,f)) and other one- 
nucleon removal nuclear processes. 
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TABLE I. Separation energies and spectroscopic fac- 
tors Sa calculated on the base of the one-body density ma- 
trix [p8|-^ for ^^O and ^^Ca. Comparison is made with 
the Hartree-Fock (HF) single-particle energies (set Ski from 
[^), natural occupation numbers Na and experimental 
data (EXP). 
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The energies are in MeV and only the states with 7 + 1 /2 are 
displayed. Experimental data are taken: "' from [BtI, ''' from 
M, '=' from il, ''' from i|. 



FIG. 1. Overlap functions (solid line), self-consistent 
Hartree-Fock single-particle wave functions (dot-dashed line) 
and natural orbitals (dashed line) for the nucleus '^"Ca. 
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